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ABSTRACT 


This  report  analyzes  and  compares  the  various  simulation  designs  and 
math  models  of  underwater  acoustic  simulation  presently  being  utilized 
in  current  training  devices  developed  by  the  Naval  Training  Device 
Center.  The  analysis  contains  model  simplifications  and  their  affects 
on  acoustic  propagation  loss.  The  comparative  analysis  includes  Device 
14E19,  Basic  AN/SQS-26  Sonar  Operator  Trainer;  Device  14B44,  P-3C  DIFAR 
Operator  Trainer;  Device  2F69,  P-3A  Weapon  System  Trainer;  Device  21A37, 
Submarine  Fleet  Ballistic  Missile  Training  Facility;  and  Device  21A39/2, 
Sonar  Room  Tactical  Team  Trainer.  The  results  of  these  models  are  then 
compared  with  the  results  of  the  ray  trace  model  developed  by  the  Naval 
Air  Development  Center. 
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SECTION  I 

ANALYSIS  OF  UNDERWATER  ACOUSTIC 

PROPAGATION  LOSS  MATH  MODELS  IN 

current  Training  devices 


K  INTRODUCTION  \ 

This  report  analyzes  and  compares  the  various  simulation  designs 
and  math  models  of  underwater  acoustic  simulation  presently  being  utilized 
in  current  training  devices  developed  by  the  NAVTRADEVCEN. 

First,  the  basic  underwater  acoustic  terms  are  defined  which  are 
required  in  the  discussion  of  the  various  math  models.  Secondly,  a 
study  and  comparative  analysis  is  made  of  the  propagation  loss  models 
utilized  in  five  current  training  devices. 

I.l  SOUND  PROPAGATION  LOSSES  Tn  THE  OCEAN^  \ 

As  sound  travels  through  the  ocean,  the  pressure  associated  with 
the  wave  front  diminishes,  this  decrease  in  pressure  is  referred  to 
as  propagation  loss.  There  are  three  factors  which  contribute  to 
this  loss. 

a.  Spreading  Loss.  Spreading  of  the  wave  front  causes  the  total 
energy  associated  with  the  wave  front  to  be  distributed  over  a  larger 
area,  resulting  in  a  decrease  in  intensity.  This  loss  is  generated 
by  the  inverse  square  spreading  law.  As  the  energy  travels  away  from 
the  source,  it  spreads  in  the  form  of  spherical  shell.  The  decrease 
in  intensity  is  exactly  proportional  to  the  increase  in  the  surface 
area  of  the  sphere.  Since  the  surface  area  is  given  by  A  =  AirR^  the 
decrease  in  intensity  is  proportional  to  the  square  of  the  radius.  The 
loss  in  dB  due  to  spreading  between  a  point  one  yard  from  the  source  and 
the  receiver  is  given  by: 

Spherical  Spreading  Loss  (db)  =  10  log  R^  =  20  log  R 

where  R  is  the  range  in  yards  between  the  source  and  receiver.  The  intensity 
loss  from  spherical  spreading  increases  about  6  dB  per  distance  doubled. 

In  cylindrical  spreading,  the  wave  front  expands  in  the  form  of  a 
cylinder  having  a  constant  height  which  is  determined  by  the  thickness 
of  the  sound  channel  or  sound  duct.  Cylindrical  spreading  loss  is  one- 
half  of  the  spherical  spreading  loss  in  dB  over  equal  ranges,  and  is 
presented  by  10  log  R  in  place  of  20  log  R.  Therefore,  intensity  loss 
from  cylindrical  spreading  increases  about  3  dB  per  distance  doubled  and 
is  given  by: 

Cylindrical  Spreading  Loss  (dB)  =  10  log  R. 
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b.  Attenuation  Loss.  The  reduction  in  sound  pressure  level  due  to 
absorption  and  scattering  is  usually  termed  attenuation  loss.  Absorption 
is  essentially  the  conversion  of  acoustic  energy  to  heat.  As  a  sound 
wave  passes  through  the  ocean  it  is  accompanied  by  successive  compres¬ 
sions  and  expansions  of  the  medium  which  produces  some  frictional  heat 
loss  from  the  sound  wave  as  well  as  heat  loss  produced  by  the  disassocia- 
tion  of  dissolved  salts  in  sea  water. 

Although  scattering  is  a  component  in  the  attenuation  of  sound,  its 
contribution  is  not  as  important  as  that  of  absorption.  When  sound  rays 
strike  reflectors  such  as  bubbles,  fish,  and  suspended  matter,  sound 
energy  is  reflected  away  from  the  direction  in  which  the  major  portion 
of  the  sound  field  is  traveling,  so  that  the  wave  itself  suffers  a  loss 
in  e'nergy  and  hence,  the  intensity  decreases. 

c.  Reflection.  The  surface  of  the  ocean  is  rarely  smooth;  therefore, 
sound  energy  striking  it  is  seldom  reflected  specularly  (mirror  reflection). 
Since  the  ocean  surface  is  constantly  changing,  the  sound  energy  is  reflected 
in  many  directions.  Surface  reflection  is  a  function  of  both  sea  state  and 
frequency. 

The  ocean  bottom  may  also  reflect  sound  waves.  Sound  reflected  from 
the  ocean  bottom  usually  suffers  a  significant  loss  in  intensity.  Part 
of  this  loss  is  due  to  scattering  effects  discussed  above.  The  ocean 
bottom  may  be  irregular  and  consequently  reflect  sound  in  many  directions. 

The  amount  of  energy  lost  through  scattering  will  vary  with  the  roughness 
of  the  bottom  and  frequency.  Most  of  the  bottom  loss,  however,  results 
from  the  fact  that  a  portion  of  the  sound  energy  will  enter  the  bottom 
and  travel  in  it  as  a  new  sound  wave.  The  amount  of  energy  which  is  lost 
into  the  bottom  will  vary  with  the  composition  of  the  bottom,  frequency 
of  the  sound  wave,  and  the  angle  at  which  the  sound  wave  strikes  the 
bottom. 

The  total  propagation  loss  over  any  given  path  is  the  sum  of  the 
losses  due  to  spreading,  attenuation,  and  reflection.  Each  of  these 
component  losses  is  subject  to  some  degree  of  variability  and  uncert- 
tainty.  For  these  reasons  the  total  propagation  loss  of  sound  waves 
traveling  between  any  two  points  in  the  ocean  cannot  be  predicted  exact¬ 
ly.  In  practice,  the  propagation  loss  over  any  particular  path  can  only 
be  estimated  by  use  of  average  measured  values  obtained  from  many  observa¬ 
tions  . 


1.2  DEVICE  I 4E I  9  PROPAGATION  LOSS  MODELS 

This  section  provides  the  propagation  loss  models  for  Basic  AN/SQS- 
26  Sonar  Operator  Trainer,  Device  14E19  (reference  1).  The  propagation 
loss  models  used  to  simulate  AN/SQS-26  surface  duct  mode  of  propagation 
are  based  upon  the  AMOS  (Acoustic  Meterological ,  and  Oceanographic 
Survey)  equations  (references  2  and  3).  AMOS  variables  and  definitions 
are  as  follows: 
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L  =  Surface  isothermal  layer  depth  (ft) 

S  =  Sea  state 

T  =  Water  temperature  (°F) 

R  =  Range  (kyds) 

Zq  *  Own  ship  depth  (ft) 
z  =  Target  depth  (ft) 
f  =  Acoustic  freq.  (kHz) 

=  Transmission  loss  (db) 
a  =  Absorption  coeff.  (db/kyd) 
ag  =  Scattering  coeff  (db/kyd) 

G  =  First  depth  loss  factor  (db) 

H  =  Second  depth  loss  factor  (db) 

The  basic  AMOS  propagation  loss  equations  are  divided  into  three  primary 
zones.  Analytical  definitions  of  these  zones  and  corresponding  propagation 
loss  equations  are  as  follows: 

a.  Direct  Radiation  Zone  (0  £  r  ^  1) 

When  both  ends  of  the  transmission  path  lie  within  or  at  the  bottom 
of  the  surface  layer,  O^Zg^l,  O^z^l,  the  following  equation  applies: 

Nw  =  20  log  R  +  aR  +  G  (z-z^)  r  +  60.  (3.1) 

At  all  other  times  the  smaller  of  the  two  transmission  losses  computed 
using  equation  (3.1)  and  equation  (3.2)  is  used. 

Nw2  =  20  log  R  +  aR  +  [25  -  /|Z  -  L|  -  /|Zo  -  L|  +  25]  (f/25)  +  60  (3.2) 

The  term  within  the  brackets  is  taken  to  be  zero  where  it  is  negative. 

b.  Zone  of  First  Order  Surface  Reflection  (zq  ^  1 ,  £  r  ^  ri  +  1/2) 

Energy  has  been  reflected  at  least  once  from  the  surface,  the  following  is 
used: 

Nw3  =  20  log  R  +  aR  +  2  (r-ri)  H  (z,  Zg)  +  [1-2  (r-r^)]  G  (z-Zg)  +  60  (3.3) 
The  shadow  zone  beyond  the  limiting  ray  is  delineated  by: 

zo  <  1,  z  <  1,  r^  ^r^(ri  +  1/2) 

For  this  zone  equation  (3.3) or  equation  (3.2)  is  used,  whichever  is  smaller. 
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c.  Zone  of  the  Second  or  Higher  Order  Surface  Reflection 

The  following  equation  is  used: 

Nw4  =  10  log  R  +  (a  +  a^)  R  +  H  (z,  zq)  -a  /I  (ri  +  ^)  +  10  log 
[/T  (r^  +  h)'\  +  60  (3.4) 

The  shadow  zone  beyond  the  limiting  ray  is  delineated  by: 

zo  1 .  z  1  1 .  >"1  +  ^  1  r 

For  this  zone,  equation  (3.4)  or  (3.2)  is  used,  whichever  yields  the  smaller 
propagation  loss. 

The  required  parameters  are  as  follows: 

(1 )  Scaled  Variables 


ri  =  %[2  -  Z  -  Zq] 


Z  <  1  ,  Zq  <  1 


Z  >  1, 


^0  >  1 


(2)  Absorption  Coefficient  (db/kyd) 

2  f  2 

■  ,  _  Af  ^T  ,  Bf 

where  A  =  .651  B  =  .0269  and  f  =  1.23  x  10%  e  +  459.6)] 

For  frequencies  less  than  4  kHz 

a  =  .-i28_fi 


4 


NAVTRADEVCEN  IH-200 


(3)  Scattering  Attenuation  Coefficient  (db/kyd) 


=  9  /~l  ■  s  >  3 

(4)  First  Depth  Loss  Factor  (db)  G  as  a  function  of  z  -  z  is 
defined  as  follows: 

G  (z  -  zo)  =  .1  X  10^'^  (z  -  Zq)  (f/25j'*^^  z  -  zo  <  1 

1/3 

G  (z  -  Zo)  =  20  (f/25)  z  -  Zq  1  1 


(5)  Second  Depth  Loss  Factor  (db)  H  as  a  function  of  z^,  and  Zq 
is  defined  as  follows: 

H  (zuzo)  =  F  (z  -  Zq)  +  F  (z)  +  F  (zq) 

where  x  is  a  dummy  variable  and  F  (x)  is  defined  as  follows: 

X  1/3 

F  (x)  =  .4  X  10  (f/8)  F  >  8  kHz 

F  (x)  =  .4  X  10^  F  <  8  kHz 


The  objective  of  this  subsection  is  to  present  the  analytic  functions  that 
are  used  to  simulate  the  AN/SQS-26  bottom  bounce  mode  of  propagation.  Bottom, 
bounce  conditions  will  be  assumed  to  exist  only  for  bottom  depths  of  6,400 
feet  or  greater  and  deflection/elevation  angles  greater  than  or  equal  to  10 
-degrees.  For  all  other  conditions,  it  is  assumed  that  propagation  takes  place 
via  surface  duct  or  shallow  water  conditions.  These  assumptions  are  based 
upon  the  recommended  modes  of  operation  given  in  reference  4.  For  the 
simulation,  it  is  assumed  that  the  propagation  loss  to  the  target  from  the 
point  of  transmission  is  equal  to  the  loss  from  the  sound  source  to  the  point 
of  reception.  Based  on  this  assumption,  the  one-way  path  loss  in  db  for  the 
bottom  path  is: 


Nw  =20  log  R  +  aR  +  Ln 

B  ^  s  s  D 


where  =  Slant  range  (yds) 
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a  =  Attentuation  coeff  (db/yd) 

Lg  =  Bottom  Loss  (db) 

®TB  =  Real  bottom  grazing  angle  (degrees) 

Initially  a  function  is  required  that  determines  the  vertical  angle  ®Bi  as 
follows:  ®Bi  =  tan'^  ^^Bi  -  °T 

where  *^Bi  =  Depth  of  bottom 
=  Depth  of  target 
=  Horizontal  range 

These  variables  are  shown  in  figure  1. 

Total  slant  range  is  calculated  by  the  following  equation: 

Rg  =  ^Bi  ^  hi  -  Pf  ^  ^hi  -  h 
^  sin  Ggi  sin  egi  sin  eei 


The  isovelocity  angle  Gg.  is  converted  to  a  real  angle  equivalent  in  order 
to  provide  the  real  transmitted  beam  deviation  loss  angles  to  the  target. 
The  transmitted  deviation  loss  angle  is  calculated  as  follows: 


0Tk  •  =  --ai  + 

Tbi  1 


/Tf^4a2  (ao  -  Gg.) 


2a' 


The  a  coefficients  were  generated  by  Honeywell  and  are  included  in  table  1. 

Gjg  is  the  real  bottom  grazing  angle  and  is  given  by  the  following  equation: 

^  1 

Cos  0jg  -  c  cos  0Jg-j 


where  c  are  constants  generated  by  Honeywell  and  given  in  table  2. 
The  attenuation  coefficient  is  given  by: 
a  =  .033  (f)^^^ 
where  f  =  frequency  (kHz) 

and  bottom  loss  simulation  functions  are  analytically  described  by  the 
following: 
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Flgvire  1.  Ray  Angles  Defined 
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TABLE  ,1.  ISOVELOCITY/REAL- RAY  ANGLE  COEFFICIENTS 


Isovelocity /Real-Rny  Angle  Coefficients 

Profile 

Depth 

(Feet) 

"o 

■'^2 

Standard 

Deviation 

0  * 
^MIN 

Degrees 

Standard 

6400 

8. 6023 

0;6603 

0.004213 

:f0.  2688 

2 

50  Foot 

9600 

9. 1153 

0. 6226 

0.004580 

iO.3446 

2 

Layer 

12,800 

8. 2957 

0. 6442 

0.004413 

.10.3275 

2 

16,000 

5.0734 

0.8200 

0.001993 

i:0.0762 

6 

25, 600 

-4. 1881 

1. 2230 

-0.002824 

10.2205 

16 

Standard 

6400 

8. 5830 

0.6575 

0.004051 

±0.3257 

2 

100  Foot 

9600 

8.8283 

0.6394 

0.004197 

±0.3230 

2 

Layer 

12,800 

8.0539 

0. 6602 

0.004199, 

±0. 3004 

2 

16,000 

4.7018 

0.8417 

0.001642 

10.0735 

6 

25,600 

-5.3533 

1. 2846 

-0.003624 

±0.3230 

16  1 
_ j 

Standard 

6400 

8. 2015 

0.6800 

0.003792 

±0.3237 

2  ! 

150  Foot 

9600 

8.5083 

0. 6594 

0.004059 

±0. 2962 

2  ! 

Layer 

12,  800 

7. 7826 

0. 6703 

0.004105 

±0. 2863 

2 

16,000 

4.3802 

0.8616 

0.001356 

±0. 2035 

6  1 

25,600 

-6.0692 

1.3195 

-0.004083 

±0.3894 

16  1 

Standard 

6400 

7.8392 

0.7091 

0.003424 

±0. 2613 

1 

4  i 

200  Foot 

9600’ 

8.0729 

0. 6832 

0.003897 

±0.2645 

4  j 

Laj'er 

12,800 

7.1932 

0.7087 

0.003510 

±0. 2254 

4  ’ 

16,000 

4. 2367 

0.8610 

0.001355 

±0.0644 

8  : 

25,600 

-3.8743 

1. 1904 

.  -0.002265 

±0. 1200 

18  ! 

Standard 

6400 

7.4898 

0. 7121 

0. 003468 

±0. 2399 

4  1 

300  Foot 

9600  ■ . 

7; 6342 

0. 6993 

0.003571 

±0.2460 

4  1 

Layer 

12,800 

6. 6991 

0.7257 

0.003225 

±0. 2212 

4  1 

16,000 

3.8759 

0.8760 

0.001217 

±0. 1382 

8  1 

25, 600 

-4.3241 

1.2033 

-0.00238,7 

±0.1158 

18  i 

1 

Standard 

6400 

6.7703 

0. 7554 

0.002769 

±0. 2747 

1 

4  ! 

400  Foot 

9600 

7.  2025’ 

0. 7099 

0.003458 

±0.2188 

1 

4 

Layer 

12,800 

6. 1839 

0.  7527 

0. 003021 

±-0.  2020 

4 

16,000 

2.  7324  ' 

0.9477 

0. 000233 

±0.3097 

8 

25, 600 

-4.5184 

1. 2034 

-0.002303 

±0. 1215 

18 

Temperate 

6400 

10,8105 

0.5667 

0.005235 

±0.3742 

2 

Summer 

9600 

10. 6152 

0.5615 

0.005308 

±0.3730 

2 

90  Foot 

12,800 

9. 6249 

0. 5991 

0.004836 

±0.3917 

2 

Layer 

16,000 

6.7419 

0. 7527 

0. 002707 

±0.1513 

4 

. 

,  25, 600 

-4.4007 

1. 2526 

-0.003285  • 

±0. 2847 

16 

Temperate 

6400 

8.4710 

0. 6676 

0.004021 

±0.2818 

'  4 

Autumn 

9600 

8. 1893 

0, 6692 

0.003965 

±0.2718 

4 

260  Foot 

12,800 

6.7723 

0.7266 

0.003110 

±0. 2283 

4 

Layer 

16,000 

3. 6068 

0.8883 

0.001142 

±0.0593 

8 

25, 600 

-5.4342 

1.2542 

-0.003012 

±0. 1717 

18 

8 
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TABLE  1.  ISOVELOCITY /REAL- uAY”  aNojlE  COEFFIC  hF.NTS  (Continued) 


Isovelocity/Real-Ray  Angle  Coefficients  (continued) 

Profile 

Depth 

(Feet) 

^0 

■  ^2 

Standard 

Deviation 

0  * 
'^MIN 

(Degi’ees) 

Temper  ate 

6400 

3.8082 

0.8795 

0.001232 

HO. 1780 

4 

Winter 

9600 

4.4706 

0.8296 

0.001974 

±0.1486 

4 

350  Foot 

12,800 

1.0820 

0.8876 

0.002349 

±0.5132 

8 

Layer 

16,000 

0.5318 

1.0083 

-0.002682 

±0. 1479 

10 

25,600 

-9.3678 

1.4144 

-0.004823 

±0.4047 

18 

Equator 

6400 

10.1425 

0. 6102 

0.004671 

±0.2882 

2 

Summer 

9600 

10. 1572 

0. 5962 

0.004847 

±0.3291 

2 

215  Foot 

12,800 

9. 0422 

0. 6379 

0.004252 

±0.3152 

2 

Layer 

16,000 

6.6005 

0.7588 

0.002645 

±0. 1591 

4 

25, 600 

-3.4573 

1. 2037 

-0.002655 

±0. 2279 

16 

Equator 

6400. 

9. 5204 

0.6509 

0.003782 

±0.2656 

4 

Autumn 

9600 

9.6214 

0. 6333 

0.004030 

±0. 2537 

4 

215  Foot 

12,800 

8.7852 

0. 6591 

0.004052 

±0.  2521 

4 

Layer 

16,000 

6. 7517 

0.7503 

0.002748 

±0. 1578 

4 

25, 600 

-2.9756 

1. 1783 

-0. 002332 

±0. 1924 

16 

Equator 

6400 

9.4627 

0. 6479 

0.004053 

±0. 2129 

2 

Winter 

9600 

9.7181 

0. 6155 

0.004542 

±0.  2463 

2 

246  Foot 

12,800 

8.9231 

0. 6462 

0.004080 

±0.3182 

3 

Layer 

16,000 

6. 6245 

0.7504 

0.002805 

±0.1439 

4 

25,600 

-3.8493 

1.2231 

-0.002927 

±0. 2506 

16 

Polar 

6400 

.  -9,.  6344 

1.4214 

-0.004949 

±0. 2494 

14 

Summer 

9600 

-13. 6549 

1.5721 

-0. 006375 

±0.4530 

16 

Zero 

12,800 

-18.4531 

1.7780 

-0.008931 

±0.5219 

18 

Layer 

16,000 

-21.4276 

1.8700 

-0.009692 

±0.4870 

20 

25, 600 

-21,9089 

1.9727 

-0.009836 

±0. 1871 

26 

Polar 

6400  • 

3.3871 

0.8592 

0.001752 

±0.1417 

4 

Autumn 

9600 

-0.3767 

1.0361 

-0. 0004649 

±0. 1561 

10 

Zero 

12, 800 

-3.6426 

1.1756 

-0.002106 

±0. 1475 

14 

Layer 

16,  000 

-6.7215 

1.2971 

-0.003489 

±0.2888 

16 

25  600 

-13.3254 

1.4774 

-0.004728 

±0.1279 

24 

Polar 

6400 

-1.9621 

1.0925 

-0.009983 

±0. 1815 

10 

Winter 

9600. 

-4. 1731 

•  1.1687 

-0.001777 

±0.1272 

'  14 

^162  Foot 

12,800. 

.  -7.4945. 

1.3072 

•  -0.003417 

±0.2077 

16  . 

Layer 

16, 000 

-10'.  7718 

1.4317 

-0.004738 

±6.  ■2356- 

18 

25, 600 

-18.8003 

1.6906 

-0.007045 

±0.1744  ^ 

24 

*  ®  minimum  ray  angle  that  will  strike  the  bottom.  Angles  less  than  0 

are  refracted  upward  before  striking  the  bottom  and  do  not  give  rise  to  bottom  bounee 
eoverage.  For  ray  angles  greater  than  50  degrees,  the  isoveloeity  and  real-ray  angles 
can  be  assumed  equal. 
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TABUS  2.  BOTTOM  BOUNCE  CONSTANTS 


• 

Bottom  Bounce  Constant  C 

Profile 

Bottom  Deptli  (Feet) 

6400 

'  9600 

12,800 

16,000 

.25,600 

STANDARD 

* 

50  Foot  Layer 

0,97504 

0.98297 

0.99346 

1.00430 

1.03882 

100  Foot  Layer 

0.97562 

0.98355 

0.99405 

1.00495 

1.03944 

150  Foot  Layer 

0.97620 

0.98114 

0.99465 

1.00555 

1.04006 

200  Foot  Layer 

0.97679 

0.98472 

0.99524 

1. 00615 

1.04067 

300  Foot  Layer 

0.97795 

0.98590 

0.99642 

1.00735 

1.04191 

400  Foot  Layer 

0.97912 

0.98707 

0.99761 

1.00855 

1.04310 

EQUATOR 

Autumn 

0.97232 

0.98082 

0.99011 

1.00099 

1.03097 

Summer 

0.97232 

0. 98102 

0,99071 

1.00118 

1.03718 

Winter 

0.97269 

0.98100 

O! 99070 

1.00119 

1.03800 

TEMPERATE 

Autumn 

0.97790 

0.98765 

0.99701 

1.00796 

1 . 04440 

Summer 

0.97268 

0.98099 

0.98931 

1,00218 

1.03841 

Winter 

0.98438 

0.99259 

1.00280 

1.01422 

1.05026 

POLAR 

Autumn 

1.00103 

1. 01301 

1.02478 

1.03613 

1,07062 

Summer 

1.02734 

1.03963 

1.05128 

1.06336  . 

1.09875 

Winter 

1,01104 

1.02250 

1.03417 

1. 04563 

1.08043 

( 
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(a)  Mud  Bottom 

l-B  *  0  ®TB  -  ^ 

Lg  “1.1  ® jg  ”  .403  ^ 

Lg  -  21 .5  ®TB  ^ 

(b)  Mud/Sand  Bottom 

Lg  =  0.1  Sjg  £  3 

Lg  =  .825  Sjg  -  2.2549  3<ejg  _<  25 


Lg  =  17.4 


®TB  ^ 


(c)  Sand  Bottom 
Lg  =  0 


®TB  -  ^ 


Lg  =  .365  Qyg  -  1.738 


5<0^g  £  36 


Lb=11.5 

(d)  Rock  Bottom 


Lg  =  .247  Sjg  -5.278 
Lb  =  7.1 

The  convergence  zone  mode  of  operation  of  the  AN/SQS-26  is  used  to  provide 
extremely  long  range  detection  capability  by  taking  advantage  of  the  refrac¬ 
tive  properties  of  the  deep  ocean.  To  determine  the  convergency  zone  range 
window  and  propagation  losses  within  the  window,  acoustic  ray  path  and  propaga¬ 
tion  loss  calculations  were  performed  and  provided  by  the  USNUSL  (Navy  Under¬ 
water  Sound  Laboratory),  (references  5  and  6).  This  data  has  been  compiled 
in  look-up  tables  and  is  not  included  in  this  report  due  to  its  voluminous 
content. 


©yg  >  36 

e  <21 
TB  - 

21<e.|.g  £  50 

"tb  "  "" 
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1.3  DEVICE  14B44  PROPAGATION  LOSS  MODELS 


This  section  provides  the  propagation  loss  models  for  the  P-3C  DIFAR 
Operator  Trainer,  Device  14B44  (reference  7)’.  The  math  model  for  total 
propagation  loss  is  given  by  the  following  equation: 

JLj  =  TL^  +  TL2  +  TL3  -TG 


where  TLt  ”  Total  propagation  loss 

TL-]  =  Loss  due  to  spreading  and  attenuation 
TL2  =  Surface  bounce  loss 
TL3  =  Bottom  bounce  loss 
TG  =  Convergence  zone  gain 

The  propagation  loss  is  determined  by  linear  approximations  of  trans¬ 
mission  loss  curves  as  shown  in  figure  2. 


In  the  surface  bounce  mode  of  propagation,  the  propagation  loss  is 
determined  as  a  function  of  sea  state  as  indicated  in  table  3.  The 
values  in  the  table  were  determined  by  extrapolating  Schulkins  empirical 
curve  (reference  8)  which  establishes  the  leakage  loss  as  a  function  of 
frequency  -  waveheight  product. 


When  operating  in  the  bottom  bounce  mode  of  propagation,  the  loss 
due  to  reflection  from  the  bottom  is  a  function  of  bottom  type  only. 

In  determining  the  bottom  loss  for  simulation,  average  reflection  values 
were  taken  for  each  of  the  three  bottom  types  simulated.  (References 
9  and  10.)  The  bottom  bounce  loss  is  assumed  to  be  zero  for  a  hard 
bottom.  For  a  medium  and  soft  bottom  the  loss  is  assumed  to  be  3  dB 
and  9  dB  respectively,  as  shown  in  table  4^ 


If  a  convergence  zone  exists  between  the  source  and  receiver,  a 
gain  factor  is  subtracted  from  the  total  transmission  loss.  The  con¬ 
vergence  zone  gain  provides  a  minimum  gain  of  12  db  with  a  further  in¬ 
crease  up  to  3  db  for  the  first  convergence  zone  and  up  to  1.5  db  for 
the  second  and  third  convergence  zones  as  a  function  of  the  distance  of 
the  hydrophone  to  the  center  of  the  zone.  The  following  model  is  based 
on  a  quantitative  discussion  of  convergence  zones  by  R.  J.  Urick  (ref¬ 
erence  9).  f.  , 

(1  -  -  Xn) 

N  ^ 


TG  =  12  + 


where  Xj  =  Target  to  sonobuoy  flat  range 
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Source  Depth:  200  Ft.  Bottom  Type: 

RCVR  Depth  :  300  Ft.  Sea  State  : 
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Figure  2.  Propagation  Loss  Device  14B44 
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TABLE  3.  SURFACE  BOUNCE  LOSS 


SURFACE  BOUNCE  LOSS  kS  A  FUNCTION  OF  SEA  STATE 

Sea  State  Sur 

face  Bounce  Loss  (TL2) 

0 

2  db 

1 

2  db 

2 

3  db 

3 

4  db 

4 

6  db 

TABLE  4.  BOTTOM  BOUNCE  LOSS 


BOTTOM  BOUNCE  LOSS  AS  A  FUNcTlON  OF  BoTtOM  TYPE 

Bottom  Type 

Bottom  Bounce  Loss  (TL3) 

Hard 

0  db 

Medi urn 

3  db 

Soft 

9  db 
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Xq  =  Ranges  between  midpoint  of  convergence  zone  and  target 
X/^  =  1/2  convergence  zone  width  • 

N  =  Number  of  the  convergence  zone 

It  is  important  to  note  that  these  models  are  limited  to  frequencies 
less  than  approximately  2500  Hz. 

1.4.  DEVICE  2F69,  PROPAGATION  LOSS  MODELS 

The  objective  of  this  section  is  to  present  the  propagation  loss 
math  model  for  Device  2F69,  P-3A  Weapon  System  Trainer  (reference  11). 

The  sound  energy  is  attenuated  in  the  simulated  environment  by  the  following 
equation: 

N  =  20  log  R  +  (  .2f  +  .00015  f^  )  R  +  60  db 

where  N  =  nominal  statistical  average  propagation  loss  between 
two  points  (db) 

where  R  =  range  from  source  (Kyd) 

where  f  =  signal  frequency  (kHz) 

The  frequency  dependent  term  is  small  for  low  frequencies  and  would 
indicate  a  20  db  loss  per  range  decade  or  again,  the  spreading  loss  varies 
as  a  function  of  the  inverse  of  R^. 

1.5  DEVICE  21 A37,  PROPAGATION  LOSS  MODELS  \ 

This  section  describes  the  acoustic  propagation  loss  model  for  Device 
21A37,  Submarine  Fleet  Ballistic  Missile  Training  Facility  (reference 
12).  The  earlier  sonar  math  modes  for  Device  21A37  were  constructed 
on  the  basis  of  ray  and  wave  acoustic  theory,  using  references  13  and 
14,  and  additional  data  was  obtained  from  reference  15.  T.  Einstein 
of  USNUSL  pointed  out  that  this  procedure  was  erroneous  since  project 
AMOS  showed  that  classical  sound  theory  applied  only  for  short  ranges 
and  at  high  frequencies.  Consequently,  work  in  this  area  was  discontinued 
and  USL  Research  Report  491  (reference  16)  was  adopted  as  the  basis  for 
sonar  simulation.  Two  separate  cases  are  considered: 

a.  Shallow  target  and  shallow  sonar 

b.  Deep  target  or  deep  sonar. 
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For  a  particular  set  of  parameters,  bottom-path  detection  range  and 
direct  path  detection  range  should  be  investigated  and  the  greater  of 
the  two  chosen.  Curve  a  represents  the  bottom  path  detection  curve; 
curve  3  the  direct-path  detection  curve  for  case  a;  and  curve  y  ,  the 
direct-path  detection  curve  for  case  b.  The  curves  for  layer  depths,  0 
ft.  and  400  ft.,  were  obtained  directly  from  the  corresponding  curves  in 
reference  16.  The  curves  corresponding  to  other  values  of  layer  depth  were 
obtained  by  graphically  interpolating  between  0  ft.  and  400  ft.,  curves. 

The  curve  for  100  ft.  layer  depth  is  shown  in  figure  3.  Each  curve  is 
individually  analyzed  and  subsequently  approximated  by  linear  segments. 


1.6-  DEVICE  2IA39/2  PROPAGATION  LOSS  MODELS 


This  section  provides  the  propagation  loss  models  for  Device  21A39/2 
Sonar  Room  Tactical  Team  Trainer  (reference  17)  Transmission  effects 
of  the  direct  path  mode  of  propagation  are  simulated  as  a  function  of 
range,  frequency,  layer  depth,  target  depth  and  own  ship  depth.  The 
total  direct  path  propagation  loss  is  taken  as  the  sum  of  a  nominal  loss 
and  a  layer  depth  correction  factor.  The  nominal  propagation  loss  is 
shown  in  figure  4,  for  representative  ranges  between  1  and  100  Kyds. 

The  layer  depth  correction  factor  A  is  defined  by  the  following  equations: 

A  =  (Z/L)  R  Z  <  2.7  L 

A  =  (2.7)  R  Z  >  2.7  L 

where  Z  =  Average  of  own  ship  and  target  depth  (ft) 

L  =  Layer  depth  (ft) 

R  =  Target  range  (Kyds) 

The  minimum  value  of  Z  is  20  ft. 

Transmission  effects  of  the  bottom  bounce  mode  of  propagation  are 
simulated  as  a  function  of  range  and  frequency  as  shown  in  figure  5. 

Bottom  bounce  loss  equals  the  value  at  4  Kyds  for  ranges  less  than  4  Kyds. 

1.7  COMPARATIVE  ANALYSIS  AND  MODEL'S IMPL I F ICATION 


This  section  includes  a  discussion  on  the  comparative  analysis  of 
the  five  previously  discussed  propagation  loss  math  models  and  the  NADC 
(Naval  Air  Development  Center)  ray  trace  model.  The  models  have  been 
analyzed  on  the  basis  of  the  major  modes  of  propagation.  Figures  6,  7, 
and  8  contain  a  comparative  analysis  of  the  direct  path  propagation,  surface 
bounce  propagation,  and  bottom  bounce  propagation  respectively.  The 
propagation  loss  curves  were  constructed  from  the  various  math  models 
as  a  function  of  the  variables  listed  in  table  5.  However,  target  and 
source  depths,  layer  depth,  temperature,  sea  state,  bottom  type  and 
depth  are  not  functions  in  all  of  the  math  models. 

Analysis  of  the  direct  path  propagation  models,  figure  6,  shows  that 
the  propagation  loss  curves  for  Devices  21A39,  21A37,  and  2F69  are  approx¬ 
imately  the  same.  Further  analysis  indicated  that  the  propagation  loss  curve 
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(qp)  S801  uo-paBScdoj^i 
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Range  (K  yds) 

Figure  3.  Propagation  Loss  Device  21A37 
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Frequence  (KC) 

Figure  4.  Nominal  Direct  Path  Propagation  Loss  Device  21A39 
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Figure  5.  Bottom  Bounce  Propagation  Loss  Device  21A39 
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1  ^ 


Range  (Kyds) 

Figure  6.  Comparative  Analysis  -  Direct  Path  Propagation 
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Figure  7*  Comparative  Analysis  -  Surface  Bounce  Propagation 
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TABLE  5.  MATH  MODEL  VARIABLES 


Math  Model  Variables 

Selected  values 

Target  depth 

200  feet 

Source  depth 

20  feet 

Layer  depth 

300  feet 

Temperature 

50°  Fahrenheit 

Sea  state 

>  3 

Bottom  type 

Soft 

Bottom  depth 

12,800  feet 

Range 

0  -  35  Kyds 

Frequency 

2.5  kHz 
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for  Device  14E19  consists  of  basically  three  separate  slopes  as  a  function 
of  range.  The  slope  for  ranges  less  than  5  Kyds  compares  favorably  with  the 
other  models  discussed.  The  major  problem  arises  from  the  slope  of  the  curve 
for  ranges  between  5  and  25  kyds  and  the  slope  for  ranges  greater  than  25  kyds. 
The  goal  to  be  realistically  achieved  would  be  +5  db  of  the  NADC  ray  trace 

model.  Therefore,  tolerance  limits  have  been  set  at  +5  db  and  are  shown  by 

the  shaded  area. 

Further  analysis  of  figure  6  indicates  definite  break  points  in  the 
propagation  loss  curve  for  Device  14B44.  The  propagation  model  consists  of 
basically  four  parts:  A  basic  loss  term,  a  convergence  zone  term,  a  surface 

loss  term,  and  a  bottom  loss  term.  The  basic  loss  term  is  derived  from 

Fleet  Numericals'  ASWEPS  data  package.  ASWEPS  is  designed  to  function  in  all 
ocean  areas.  Incoming  synoptic  reports  are  summarized,  analyzed,  and  then 
transmitted  daily  from  the  U.S.  Naval  Oceanographic  Office.  Passive  sonar 
ranges  are  derived  from  a  high  speed  propagation  loss  model  valid  for  low 
frequencies.  This  model  is  a  modified  ray  trace  which  calculates  losses 
via  surface  duct,  and  convergence  zones,  assuming  initially  that  the  ocean 
is  isovelocity  and  then  making  an  empirical  correction  for  refraction. 

The  ASWEPS  data  is  broken  down  for  various  ocean  areas,  summer  and  winter 
conditions,  and  a  combination  of  four  source  and  receiver  depths.  This 
information  is  provided  at  four  discrete  frequencies  up  to  1,700  Hz.  The 
NADC  ray  trace  model  was  chosen  for  comparison  here  since  it  is  very  similar 
to  the  ray  trace  model  utilized  in  generating  the  ASWEPS  curves. 

In  the  Device  14B44  propagation  loss  model,  all  source  and  receiver  depths 
have  been  averaged,  then  the  data  is  averaged  over  winter  and  summer  and 
finally  averaged  over  all  areas.  The  result  is  four  curves,  one  for  each 
frequency,  of  propagation  loss  versus  range.  These  curves  are  then  approx¬ 
imated  with  four  linear  approximations  resulting  in  the  definite  break 
points  in  the  curves.  There  are  also  very  significant  differences  in 
propagation  loss  between  different  areas  and  also  in  the  same  area  for 
different  seasons  which  have  an  important  bearing  on  training,  as  illustrated 
by  the  difference  in  the  propagation  loss  resulting  from  the  averaging  of 
depths,  seasons,  and  latitudes  in  the  Device  14B44  model. 

Analysis  of  the  surface  bounce  propagation  models,  figure  7,  again 
shows  the  steeper  slope  of  the  Device  14E19  propagation  loss  curve  and 
the  definite  break  points  in  the  Device  14B44  propagation  loss  curve. 

Although  the  basic  loss  term  for  the  Device  14B44  propagation  loss  model 
is  derived  from  the  ASWEPS  data,  the  surface  loss  term  is  derived  from  AMOS 
data.  The  data  from  which  the  AMOS  results  were  derived  only  covered  the 
frequency  range  from  2  KHz  to  25  KHz.  The  AMOS  results  were  not  intended 
to  be  used  at  low  frequencies.  To  extend  these  results  below  2  KHz  will 
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result  in  error.  An  error  exists  in  the  fact  that  the  same  surface  bounce 
loss  is  added  to  the  basic  curves  at  all  ranges.  In  the  real  world  the 
surface  bounce  Toss  increases  as  a  function  of  range  and  the  number  of 
surface  bounces.  Still  another  error  exists  in  that  the  surface  loss  has 
previously  been  included  in  the  composite  ASWEPS  data  and  should  not  be 
included  a  second  time. 

Analysis  of  the  bottom  bounce  propagation  models,  figure  8,  reveals 
an  unusually  large  initial  propagation  loss  variation.  The  deviation  here 
can-be  as  much  as  25  dB  and  should  be  revised  by  adding  a  correction  factor 
for  short  ranges. 

The  bottom  bounce  propagation  loss  model  for  Device  14E19  was  based 
on  the  assumption  of  ideal  straight  line  (isovelocity)  propagation  paths. 

In  the  real  world,  the  ocean  rarely  provides  ideal  isovelocity  conditions 
and  the  propagation  paths  differ  radically  from  straight  line  paths.  The 
variation  depends  on  the  ocean  velocity  structure.  Even  though  real  ray 
angles  are  computed  from  the  equivalent  isovelocity  angles,  it  appears 
as  though  an  error  is  introduced  here,  particularly  at  short  ranges. 

In  the  Device  14B44  bottom  bounce  propagation  model,  if  the  bottom 
loss  is  considered  important,  a  loss  of  0,  3,  or  9  dB  is  added  to  the 
basic  loss  term  depending  on  the  bottom  type.  The  same  bottom  loss  is 
added  to  the  basic  curves  at  all  ranges.  This  is  contrary  to  what  ac¬ 
tually  happens.  Operating  in  deep  water  at  low  frequencies,  energy 
incident  on  the  bottom  at  low  grazing  angles  suffers  considerably  less 
than  higher  grazing  angles.  The  same  error  also  arises  here  as  in  the 
surface  loss  propagation  model,  in  that  the  bottom  loss  has  previously 
been  included  in  the  composite  ASWEPS  data  and  should  not  be  included 
a  second  time. 

To  aid  in  error  analysis,  relative  error  v/as  plotted  versus  range  for 
each  model  and  each  major  mode  of  propagation  as  shown  in  figures  9,  10, 
and  11.  A  positive  error  indicates  a  decreased  propagation  loss  resulting 
in  a  stronger  signal  strength  and  premature  target  detection.  A  negative 
error  indicates  too  great  a  magnitude  of  propagation  loss  resulting  in  a 
weaker  signal  and  delayed  target  detection.  At  short  ranges  it  is  shown 
that  the  error  is  consistently  positive,  resulting  in  premature  target 
detection.  At  greater  ranges,  it  is  found  that  the  relative  error  is 
still  generally  positive  although  considerably  less  than  that  found  at 
shorter  ranges. 
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Figure  9.  Relative  Error  —  Direct  Path  Propagation 
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Figure  lO,  Relative  Error  --  Surface  Bounce  Propagation 
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APPENDIX  A 
SYMBOL  DEFINITIONS 
L  =  Surface  isothermal  layer  depth  . 

S  =  Sea  state 
T  =  Water  Temperature 
Rh  =  Horizontal  Range 
Rg  =  Slant  Range 
Z  =  Ownship  depth 
z  =  Target  depth 
f  =  Acoustic  frequency 
Nw  =  Transmission  loss 
a  =  Absorption  coefficient 
as  =  Scattering  coefficient 
G  =  First  depth  loss  factor 
H  =  Second  depth  loss  factor 
r,  ri ,  z,  zo  =  Scaled  variables 

Og-j  =  Isovelocity  range  angle 
0TBi  =  Transmitted  deviation  loss  angle 
QjB  =  Real  bottom  grazing  angle 
^Bi  =  Depth  of  bottom 
=  Depth  of  target 

Lg  =  Bottom  bounce  loss  as  function  of  bottom  type 
TLj  =  Total  propagation  loss 
TLi  =  Loss  due  to  spreading  and  attenuation 
TLz  =  Surface  bounce  loss 
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APPENDIX  A  (Continued) 

TL3  =  Bottom  bounce  loss 
TG  =  Convergence  zone 

=  Target  to  sonobuoy  flat  range 
Xg  =  Range  between  midpoint  of  convergence  zone  and  target 
X/\  =  1/2  convergence  zone  width 
N  =  Number  of  convergence  zone 
A  =  Layer  depth  correction  factor 


( 
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APPENDIX  B 

FLOW  DIAGRAM  AND  FORTRAN  IV  PROGRAM 

Appendix  B  contains  the  flow  diagram  and  FORTRAN  IV  Program  for  Device 
14E19  math  model.  The  FORTRAN  Program  was  generated  by  J.  H.  Gardner 
(Code  38)  to  aid  in  calculating  the  Device  14E19  math  model  from  the 
following  input  data;  target  depth,  layer  depth,  own  ship  depth,  actual 
range,  acoustic  frequency,  and  sea  state.  With  this  program  many  sets 
of  data  can  be  computed  on  the  Sigma  7  computer  in  a  reasonably  short 
time  for  comparative  analysis  with  other  math  models. 
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CALCULATION  OF  PROPAGATION  LOSS  FOR  DEVICE  14E19 

J.  H.  GARDNER  -  CODE  381 

SURFACE  BOUNCE  WITH  SHADOW  ZONE 
DIRECT  PATH  WITH  SHADOW  ZONE 
DATA  INPUT  AND  VARIABLES  WILL  BE  READ  AS 
HORIZONTAL  RANGE  =  HRANGE 
SLANT  RANGE  =  SLANTR 
BIG  Z  =  TARGET  DEPTH  =  TDEPTH 
BIG  L  =  LAYER  DEPTH.  =  XDEPTH 
BIG  ZO  =  OWN  SHIP  DEPTH  =  ODEPTH 
‘^FA  ‘^TATF  = 

F  =  ACOUSTIC  FREQUENCY  =  ACFREQ 
BOTTOM  SLOPE  =  BSLOPE 
DEPRESSION  ANGLE  =  DEPRAG 
BOTTOM  DEPTH  =  BDEPTH 
SCALED  RANGE  =  SRANGE 
ABSORPTION  COEFFICIENT  =  A 

ACTUAL  RANGE  =  ARANGE 
SCALED  OWNSHIP  DEPTH  =  ZO 
SCALED  TARGET  DEPTH  =  Z 
SHADOW  ZONE  SCALE  FACTOR  =  SZSF 


FIRST  DEPTH  -  LOSS  FACTOR(DB)  =  G 
SECONDDEPTH  -  LOSS  FACTOR(DB)  =  H 
FORMAT  (4F6.2,  F7.2,ll) 

READ  (105,  5)  TDEPTH, XDEPTH, ODEPTH .ARANGE, ACFREQ, JSEA 
IF  (TDEPTH  .EQ.  999.99)  IF  (XDEPTH  .EQ.  999.99) 

$STOP 

COMPUTE  ABSORPTION  COEFFICIENT***COMPUTE  SCALE  OWNSHIP  DEPTH 
COMPUTE  SCALE  TARGET  DEPTH***COMPUTE  SCALE  RANGE 
A  =  0.033*ACFREQ  **(3./2) 

ZO  =  (0DEPTH/XDEPTH)**.5 
Z  =  (TDEPTH/XDEPTH)**.5 
SRANGE  =  ARANGE/(XDEPTH)**.5 

IF  (Z<1.)IF(Z0<1.)  SZSF=((l.-Z0)/4.)+(SQRT((Z**2.)-l))/5.); 
GO  TO  3 

SZSF  =  (SQRT((Z0**2.)  -  1.)  +  SQRT  ( (Z**2. )-l . ) )/5 . 

COMPUTE  FIRST  DEPTH  -  LOSS  FACTOR 
IF  (Z  -ZO  <  1.)  GO  TO  8 
G  =  ((ACFREQ/25. )**(!. /3.))*20. 

GO  TO  9 

G  -  ((.1)*(10.)**(2.3*(Z  -  Z0)))*( (ACFREQ/25.)**. 3334) 
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9  IF  (ACFREQ  <  8.)  GO  TO  10 

COMPUTE  SECOND  DEPTH  LOSS  FACTOR 

H  =  .4*((ACFREQ/8. )**(!. /3.))*((10.**(Z-Z0))  +  (10.**Z)  +(10.**Z0)) 
GO  TO  11 

10  H  =  .4*((10.**(Z-Z0))+(10.**Z)+(10.**Z0)) 

11  IF  (JSEA  <  3)  GO  TO  12 

COMPUTE  SCATTERING  ATTENUATION  COEFFICIENT 
AS  =  ((ACFREQ/XDEPTH)**.5)*9. 

GO  TO  14 

12  AS  =  ((ACFREQ/XDEPTH)**.5)*4.5 

BASIC  AMOS  PROPAGATION  LOSS***START  WITH  XNW2 

14  XYZ  =  25.-((ABS(TDEPTH-XDEPTH))**.5)  -  ((ABS(0DEPTH-XDEPTH))**.5) 
1+5.  *  ARANGE 

IF  (XYZ)15,16,16 

15  XYZ  =  0.0 

16  XNW2  =  20.  *  (ALOGIO(ARANGE))  +  (A  *  ARANGE)  +  XYZ  *( (ACFREQ/25 . ) 
l**(l./3.))  +  60. 

COMPUTE  DIRECT  RADIATION  ZONE 
IF  (SRANGE  .GE.  0.)IF(SRANGE  <  SZSF)  GO  TO  20 
IF(SZSF  <  SRANGE)  IF  (SRANGE  <  SZSF  +.5)  GO  TO  21 
IF  (SZSF  +  .5  <  SRANGE)  GO  TO  22 

20  XNWl=20.*AL0G10(ARANGE)+(A*ARANGE)+(G*(SRANGE/SZSF)+60.) 

IF  (ZO  >  0.)  IF  (Z  >  0.)  IF  (ZO  <  1.)  IF  (Z  <  1.) 

$WRITE  (108,30)XNW1 ;  GO  TO  7 
IF  (XNW2  >  XNWDWRITE  (108,30)XNW1 ;  GO  TO  7 
WRITE  (108,31 )XNW2 
GO  TO  7 

21  XNW3  =  20.*(AL0G10(ARANGE))  +  (A*ARANGE)+  2.*(SRANGE-SZSF)*H+ 
$(1.-2.*(SRANGE  -  SZSF))*  G  +  60. 

IF  (ZO  >  1.)  WRITE  (108,36);  GO  TO  7 
IF  XNW2  >  XNW3)WRITE(108,32)XNW3;  GO  TO  7 
WRITE  (108,33)XNW2 
GO  TO  7 

22  XNW4  -  10.  *  (AL0G10(ARANGE))+(A+AS)*  ARANGE  +  H 
1-AS  *  ((XDEPTH  **.5)  *  (SZSF  +  .5))  +  10.  * 

2(AL0G  10((XDEPTH  **  .5)  *  (SZSF  +  .5)))  +  60. 

IF  (ZO  >  1.)  WRITE  (108,37);  GO  TO  7 

IF  (XNW2  >  XNW4)  WRITE  (108,34)XNW4;  GO  TO  7 
WRITE  (108,35)  XNW2 
GO  TO  7 

30  FORMAT  (3/30H  WITHIN  DIRECT  RADIATION  ZONE  ,//12X,8HNW(l )  =  ,F6.2, 
$4H(DB)) 

31  FORMAT  (3/30H  WITHIN  DIRECT  RADIATION  ZONE  ,//12X,8HNW(2)  =  ,F6.2, 
$4H(DB)) 

32  FORMAT  (3/62H  WITHIN  ZONE  OF  FIRST  ORDER  SURFACE  REFLECTION  AND  SHA 
$  DOW  Z0NE//12X,8HNW(3)  =  ,F6. 2 ,4H(DB) ) 

33  FORMAT  (3/62H  WITHIN  ZONE  OF  FIRST  ORDER  SURFACE  REFLECTION  AND  SHA 
$DOW  Z0NE//12X,8HNW(2)  =  ,F6. 2 ,4H(DB) ) 

34  FORMAT  (3/  'WITHIN  ZONE  OF  SECOND  OR  HIGHER  ORDER  SURFACE  REFLECT 
$ION  Z0NE'/12X,8HNW(4)  =  ,F6. 2 ,4H(DB) ) 
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35  FORMAT  (3/  'WITHIN  ZONE  OF  SECOND  OR  HIGHER  ORDER  SURFACE  REFLECT 
$ION  Z0NE'/12X,8HNW(2)  =  ,F6.2,4H(DB) ) 

36  FORMAT  (3/  'ZO>  1, WITHIN  ZONE  OF  FIRST  ORDER  SURFACE  REFLECTOR  AND 
1  SHADOW  ZONE-NW3') 

37  FORMAT  (3/  'ZO  >  1,  WITHIN  ZONE  OF  SECOND  OR  HIGHER  ORDER  SURFACE  R 
lEFLECTION  Z0NE-NW4') 

END 
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